Spin-dependent inertial force and spin current in accelerating systems 
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The spin-dependent inertial force in an accelerating system under the presence of electromagnetic 
fields is derived from the generally covariant Dirac equation. Spin currents are evaluated by the 
force up to the lowest order of the spin-orbit coupling in both ballistic and diffusive regimes. We 
give an interpretation of the inertial effect of linear acceleration on an electron as an effective electric 
field and show that mechanical vibration in a high frequency resonator can create a spin current via 
the spin-orbit interaction augmented by the linear acceleration. 



I. INTRODUCTION 

Studies of inertial effects on electrons have a long his- 
tory dating back to 1910s. Barnett investigated the mag- 
netization induced by rotational motion. 1 Einstein and de 
Haas carried out the reversed experiment. 2 They mea- 
sured the gyromagnetic ratio and the anomalous g factor 
of electrons before the establishment of modern quan- 
tum physics. Stewart and Tolman estimated the elec- 
tron mass by measuring the charge accumulation at the 
rim of a metal due to the linearly accelerating motion. 3 
Rapid progress in nanotechnology has allowed us to study 
the coupling of mechanical motion and electromagnetism 
in the quantum mechanical regime. Effects of mechan- 
ical rotation on nanostructured magnetism are detected 
in microcantilevers 4 ' 5 and a torsional resonator 6 . The 
quantization of the rotational motion is observed in mag- 
netic nanoparticles. 7 Single quantum excitations of vibra- 
tion, namely, phonons, can be controlled in a piezoelec- 
tric acoustic wave resonator. 8 There has been theoretical 
work on effects of mechanical rotation on nanostructured 
magnetism, 9-15 the coupling of nanomechanical vibration 
and magnet, 16-19 and the phonon-spin coupling related 
to spin relaxation. 20-26 However, the contribution of the 
spin-orbit interaction (SOI) in accelerating frames has 
not thoroughly been studied in the previous papers. 

Recent developments in spintronics, 27 ' 28 which relies 
on not only electron's charge but also its spin, have en- 
abled us to utilize a "spin current" , a flow of spins. In 
this context, the coupling of magnetization and spin cur- 
rent is of great interest in the field of spintronics. 29-32 To 
harness the spin current, the understanding of the spin- 
orbit interaction is indispensable. Recently, the present 
authors proposed a fundamental theory describing the di- 
rect coupling of the mechanical rotation and spin current 
and predicted the spin current generation arising from ro- 
tational motion. 33,34 Our finding has made the first step 
to extend the theory of spin current in the inertial frame 
to that in the nonincrtial frame. In this paper, we pro- 
vide a systematic approach to study spin current genera- 
tion from mechanical motion, including time-dependent 
rigid rotation and linear acceleration. Firstly, we derive 
the spin-dependent inertial force in a rotating frame in 



the presence of an applied magnetic field. This force is 
responsible for the generation of a spin current due to 
mechanical rotation. Secondly, generation of the spin 
current by linear acceleration is investigated. We show 
that nanomechanical vibration can create an ac spin cur- 
rent on the basis of the inertial spin-orbit coupling in 
accelerating systems. 

The outline of the paper is the following. In Sec. II, 
we review the Dirac equation in a rotating frame. In 
Sec. Ill, the Pauli-Schrodinger equation in the rotating 
frame is derived. In Sec. IV, we derive a full expres- 
sion of a spin-dependent force caused by mechanical ro- 
tation. Spin current generation in the presence of impu- 
rity scattering is studied in Sec. V. A detection method 
of such spin currents is proposed in Sec. VI. Spin cur- 
rent generation from mechanical vibration is investigated 
in Sec. VII. The paper ends with a few concluding re- 
marks in Sec. VIII and three Appendices. Appendix A 
contains short summaries on vierbein. Details of Foldy- 
Wouthuysen-Tani transformation are given in Appendix 
B. Electromagnetic fields in a rotating frame are briefly 
summarized in Appendix C. The spin diffusion equation 
in a rotating frame is solved in Appendix D. 



II. DIRAC EQUATION IN A ROTATING 
FRAME 

In this section, we review the Dirac equation in a ro- 
tating frame. According to Einstein's principle of equiv- 
alence, gravitation cannot be locally distinguished from 
inertial effects due to acceleration of the frame of refer- 
ences. In the general relativity, both gravitational and 
inertial effects are expressed by a metric and connection 
in a curved space-time. The fundamental equation of a 
spin- 1/2 particle in a curved space-time is the generally 
covariant Dirac equation : 
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where c,h,q — — e, and m are the velocity of light, the 
Planck constant, the charge and mass of an electron, 
= (Aq,A) is the U(l) gauge potential, and T M the 
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spin connection 35 . The Clifford algebra in the curved 
space-time 7 M = 7^(2;) satisfies 

{r(x),7 v (x)} = 2g^(x) (2) 

where g^ v (x) (/j,,v = 0,1,2,3) is the inverse of the 
coordinate-dependent metric g^ v (x). The coordinate 
transformation from a rigidly rotating frame to the in- 
ertial frame is given by 

dr' = dr+(flx r)dt, (3) 

where the rotation frequency with respect to an inertial 
frame is fl(t). Here we assume that the rotation velocity 
ft x r is much less than the speed of light. The space-time 
line element in the rotating frame is 

ds 2 = g^dx^dx" 

= [-c 2 + (ftx r) 2 ]dt 2 + 2(ft x v)dtdr + dr 2 . 

Thus, the metric in the rotating frame becomes 

( -l + u(x) 2 u x {x) Uy(x) u z (x) 
u x {x) 10 

u y (x) 10 

u z {x) 1 



9^ 



\ 



(4) 



(5) 



with 



(6) 



u(x) = fl(t) x r/c. 

This metric leads to the Clifford algebra and the spin 
connection in a rotating frame as 

7 (x) = if3, 7*(x) = iftaii - Ui(x 

ft ■ s 

r = 



where 
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r t = 0, 



(7) 
(8) 



O a- 
er O 



(9) 



are the Dirac matrices and S is the spin operator for 
4-spinor defined by 



£ = —a x a 

4:1 



a O 
O er 



(10) 



with the Pauli matrix er (details of the spin connection 
are given in Appendix A). From Eqs. (l)-(8), the Dirac 
equation in a rotating frame is written as 

£T*, (11a) 



ih 



dt 

H = I3mc 2 + ca ■ n + qA Q - ft ■ (r x tt + S) , (lib) 

where n = p — qA is the mechanical momentum and r is 
the position vector from the origin at the rotation axis. 

In classical mechanics, the Hamiltonian in the rotating 
frame has the additional term fl ■ (r x tt) which repro- 
duces the inertial effects: Coriolis, centrifugal, and Eu- 
ler forces 36 . The term ft ■ S is called the spin- rotation 
coupling. 37 " 39 The last term of Eq. (lib), ft ■ (r x tt + S), 
can be regarded as a quantum mechanical generalization 
of the inertial effects obtained by replacing the mechani- 
cal angular momentum r x tt with the total angular mo- 
mentum r x tt + S . 



III. PAULI— SCHRODINGER EQUATION IN A 
ROTATING FRAME 



The Dirac equation is an equation of 4-spinor wave 
function, which contains the up/down-spin electron and 
positron components. As the energy gap between the 
electron and positron state is much larger than the en- 
ergy level of condensed matter systems, we take the 
low energy limit of the Dirac equation to obtain the 
Pauli-Schrodinger equation of up- and down-spin elec- 
trons. Following the low energy expansion and block- 
diagonalization method of the Dirac Hamiltonian devel- 
oped by Foldy, Wouthuysen, 40 and Tani, 41 we derive the 
Pauli-Schrodinger equation in a rotating frame (see Ap- 
pendix B for the detail of the derivation). The Hamil- 
tonian (lib) is divided into the block diagonal and off- 
diagonal parts denoted by £ and O, respectively; 



with 



H = f3mc 2 +£ + (D, (12) 

£ = qA -ft-(r xtt + S), (13a) 
O = ca ■ tt. (13b) 

By successive Foldy- Wouthuysen- Tani transformations, 
the Hamiltonian up to the order of 1 /m 2 becomes 



H=f3 



mc + 



O 2 



2mc 2 



+ £ 



8to 2 c 4 



O, {€>,£} + ihO 



(14) 

Neglecting the rest energy term in Eq. (14), the Pauli- 
Schrodinger equation for the upper component of Dirac 
spinors, namely, 2-component electron wave function, in 
the rotating frame is obtained as 33 



where 



with 



and 



i?l ~dt = HvK ^'' 

Hpu = H K + H z + Hi + H s + H D , 



H K = ^-TT 2 + qA , 

Hz = fJ-BC ■ B, 

Hi = -ft ■ (r x tt + S), 

H s = |^<T • (tt x E' - E' x tt), 
H D = -^divE', 



i =— A = - S = -<x 

2m ' 4m 2 c 2 2 



(15) 
(16) 

(17a) 

(17b) 
(17c) 

(17d) 
(17e) 



(18) 



E' = E + (O x r) x B. 



(19) 

The Hamiltonian in Eq. (15), i/pR, is a 2 x 2 matrix 
operator and ip is the 2-spinor wave function of a single 
electron. 
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A. Lowest order of the expansion 

In the lowest order of the expansion, the Hamiltonian 
to the order of 1/m is given by H^ + Hz + Hi. The 
spin- independent H^ contains the kinetic energy and the 
potential energy. The Zeeman energy Hz contains the g 
factor of the electron equal to 2. Combined H^ with Hz, 
the coupling with magnetic field, 



^(rx 7r + 2S) 



B 



(20) 



is obtained, which contrasts with Eq. (17c): the mechan- 
ical rotation couples to the total angular momentum of 
the electron 



r x 7r + S. 



(21) 



The inertial effects, namely, the Coriolis, centrifugal, and 
Euler forces, are reproduced by the first term of Hj as 
mentioned above. The second term of Hi is the spin 
rotation coupling term. Introducing the "Barnett field" , 



Bo = (m/g)0, 



(22) 



we can combine the spin rotation coupling with the Zee- 
man term, leading to a different form: 



H B tr -(B + B a ), 



(23) 



Equation (23) shows that the spin-rotation coupling can 
be interpreted as a correction to the Zeeman effect with 
an effective magnetic field Bn- 

Previous theoretical work 9-15 has been done on the 
basis of the Barnett field. In the following sections, we 
study inertial effects on spin current using the spin-orbit 
interaction H$ which is obtained from the second order 
of the expansion. 



B. Second order of the expansion 

The expansion of the order of 1/m 2 yields the SOI 
and Darwin term with the mechanical rotation, Hg and 
Hi). In the absence of the rotation, ft = 0, these terms 
reproduce the conventional SOI and Darwin terms in the 
rest frame. In the presence of the rotation, we find that 
the electric field E in the two terms in the inertial frame is 
modified by an additional term (17 x r) x B. This result 
is consistent with a general coordinate transformation 
of electromagnetic fields between the rest frame and the 
rotating frame 42 (the derivation of the transformation is 
given in Appendix C). 



C. Renormalization of SOI 

The contribution of H$ to Hi in vacuum is negligible 
since the dimensionless spin-orbit coupling parameter 



VSO = 



X(mv) 2 



= (- 



<C 1 



(24) 



with the electron velocity v. However, the spin-orbit cou- 
pling is enhanced in metals and semiconductors such as 
p^. 45,46 rp ne renorma ]j za tj on depends on detailed elec- 
tronic structures and electron correlations. 47,48 In the 
present study, we do not go into detail about this proce- 
dure. Nevertheless, the results obtained in this paper are 
universal in nature and one can start with an effective 
Hamiltonian, e.g., Luttinger 49 or Rashba 50 model, which 
refers the electromagnetic fields (19) in a rotating frame. 
Replacing the momentum mv with the Fermi momentum 
Hkp, the coupling 77 SO becomes \k 2 F where A is an en- 
hanced spin-orbit coupling parameter. The coupling 7730 
of Pt is estimated as 0.59 by the nonlocal measurement of 
the spin Hall effect. 45 ' 46 Electrons in a noninertial frame 
cannot distinguish the inertial effect originating from the 
mechanical rotation (ft x r) x B in Eq. (17d) from con- 
ventional electric field E. Thus, the new effect due to the 
SOI with the mechanical rotation can be sizable effects in 
the large SOI systems as shown in the following sections. 



IV. SPIN-DEPENDENT INERTIAL FORCE IN 
A ROTATING FRAME 

Let us consider semi-classical equations of motion for 
an electron based on the Pauli-Schrodinger equation in 
a rotating frame. A quantum mechanical analogue of a 
"force" T is defined by 



1 <9r 
jr = _ [mi ., HpR }+m-, 



(25) 



with r = [r, Hpn]/ih. FromEq. (15), the spin-dependent 
velocity including the effect of mechanical rotation is ob- 
tained as 



with 



V + V! + V CT , 



v=l[r,ff K ] = -, 

in m 

vi = — [r,Hi] = ft x r, 
in 

v CT = \[r,H s ] = ^-ctx E'. 
in n 



(26) 

(27a) 
(27b) 
(27c) 



The "force" T = + T\ + F 2 + T t is obtained as 

Jo =q[E' + v x (B + 2B n )] + mft x (ft x r), (28a) 

Ti =-^-{(<x x E') x (B + Bo) - [(B + B n ) x a] x E'} 



qm\ 



a ■ (ft x v)B + 2(B • fl)cr x v - (B • v)cr x £1 



+n ■ (r x B)cr x O — (B • fi)er x (ft x r)], 



(28b) 



J 2 = 



mq 2 \ 2 



h 



(a ■ E')mv x E' 



+i[((T ■ E')B x ft - (<t ■ B)E' x ft + (B • ft)E' x tr] 



+(E' x B) x fl + 2(B-ft)E' 
dft qmX 



F t =mr x — 
at 



-a x 



r x — — x B 

dt 



(28c) 
(28d) 



4 



The spin-independent To consists of the electromagnetic 
force in a rotating frame q(E' + v x B), Colioris force 
qv x 2Bq = 2mv x and centrifugal force mfl x (f2 x r). 
The first term in T% is the Euler force. The other terms 
in T with the spin operator er are responsible for spin- 
dependent transport of electrons. The full expression of 
the spin-dependent force in a rotating frame in the pres- 
ence of electromagnetic fields is one of the principal re- 
sults of this paper. In the absence of rotation, Q = 0, the 
above expression of T reproduces the previous results in 
an inertial frame. 51 

In order to have the better understanding of the spin- 
dependent force, we study the case of = (0, 0, 0), B = 
(0,0,S),E = 0, |B/Bn| > 1 and neglect the terms of 
the order of 7yg , |fi/aj c | 2 with the cyclotron frequency 



qB 



(29) 



T is decomposed into the xy- and z-components, T x and 
Tn . Thus, we have 



with 



T x ~ q(K + E CT + v x B) 



E r = (B • n)r, 

E CT = -^(B-o-)(B-fi)r. 



(30) 

(31a) 
(31b) 



Here, E r is an electric field induced in the rotating frame 
with an applied magnetic field B, and E CT is an "effective 
spin-dependent electric field" induced by the SOI, Hs- 
Coexistence of electric and magnetic fields, E and B, 
yields the E x B drift, which is the motion of guiding 
center of a charged particle with the drift velocity 43 



v d = 



E x B 
~ B2 - ' 



(32) 



From Eq. (30) we obtain two types of the drift motion for 
an electron wave packet: one is the charge drift motion 
with 



E r x B 

B 2 ' 



(33) 



and the other is the spin-dependent drift motion with 

E„xB 



vt = 



B 2 



(34) 



Figure 1 (a) illustrates the relation of the rotation, mag- 
netic field, induced spin-dependent field, and drift veloc- 
ity. In a ballistic regime, the latter produces the spin 
current in the azimuthal direction, 



J, = en.Trcr,v„ = 2nenuj r Rej 



(35) 



where R is the distance from the rotation axis, the 
unit azimuthal vector, n the electron density, and the 
dimensionless parameter 



Afcp 



(36) 





FIG. 1. (a) Spin-dependent electric field E ff and drift ve- 
locity v^. are illustrated. An external magnetic field B is ap- 
plied along the rotation axis (z-direction). For the z-polarized 
spins, the electric field, E+(E_), is induced in the radial out- 
ward(inward) direction, (b) The drift velocities in oppo- 
site directions result in the spin current jf in the azimuthal 
direction. 



with Fermi energy cf 33 Setting B = IT, £1 = lkH, 
Xkp ps 0.6, kp ~ 10 10 m, and R = 10mm, |J S | is esti- 
mated to be about 10 s A/m 2 . 

Next, we consider fluctuation of the rotating axis, 
caused by a high speed rotor in operation. Such fluc- 
tuation effects on spin current come from f2 dependence 
of E CT in Eq. (31b) as well as the time-derivative of the 
f2 in the second term of Tt ■ If the rotation frequency has 
a time-dependent component as 



where 



n(t) = n + sn{t) 



sn(t) = 6to x (t) + 6n«(t), 



(37) 



(38) 



with Sftx -L Oo, 5Sl\\ || fi , and \Sfl\ < |0| , time- 
dependent drift motion is induced: 



<5E CT x B 

B 2 



with 



SE^ = -^(<7-B)(B-<yn)r, 

n 



(39) 



(40) 



and the fluctuation of the spin current is obtaind as 

r . , 2ne 2 XBR r , , 

63.(t) = Ml||(t)e*. (41) 

The time-derivative of Eq. (38) is divided into the xy- 
and z-components: 



dsn _ /asn\ /dsn\ 

~df ~ \~df) ± + \~dfju ' 



(42) 



From the second term of Tt , we have an additional effec- 
tive spin-dependent electric field created in the xy-plane 
SE' as 



qm\ 



fdSfl\ 

B)rx hr 



(43) 



This yields spin current in the azimuthal direction 

2nem\R ( dSfl\ 
5J' s (t) = — E — ST e ^ ( 44 ) 



The ratio \8J'J5J S \ is equal to | ^ | /\5Clw c \. The time 
scale of the fluctuation of rotation is usually much smaller 
than that of the cyclotron frequency. Thus, the contribu- 
tion from Tt to the spin current fluctuation is negligible. 




V. EFFECTS OF IMPURITY SCATTERING 

In this section, we discuss the spin current generation 
in the presence of impurity scattering. We consider a 
Pt thin film attached on a rotating disk with a uniform 
magnetic field parallel with the rotation axis as shown 
in Fig. 2. Spin dependent transport in a system with a 
strong SOI can be described by semiclassical equations 52 



r = v + v ct0 , 
hk = -e (E + f x B) . 



(45a) 
(45b) 



Here, v = Tik/m represents the normal velocity, k is the 
wave vector, E and B are applied electric and magnetic 
fields. The anomalous velocity ivo originating from the 
SOI in an inertial frame is written as 



eX 



v ct0 



-<x x E. 



(46) 



It is straightforward to extend the equations in an in- 
ertial frame to that in a rotating frame using Eqs. (25) 
and (26), 



Vi 



Kk = T. 



(47a) 
(47b) 



Let us consider spin current generation in a rotating nor- 
mal metal with the large spin-orbit coupling such as Pt 
in the presence of spin-independent impurity scattering. 
The electron distribution function depends on spins be- 
cause of the spin-dependence of the semi-classical equa- 
tions (47). In this case, the transport equation of non- 
equilibrium steady states is written as 



dv 



<9k 



fcr — /o 



(48) 



where / CT = J CT (r,k) is the spin-dependent distribution 
function, /o = /o(e) the Fermi-Dirac distribution func- 
tion, and r the relaxation time. 

Combining the semiclassical equations (47) and with 
Eq. (48), by putting E = and 



T = T \ 



-e[E r 



E„ 



xB], 



the solution is 



U = /o + ev • r 



E; + ro; e x E^ df 
1 + (rw c ) 2 de '■ 



(49) 



(50) 



FIG. 2. The z-polarized radial spin current is converted to 
the inverse spin Hall voltage in the azimuthal direction e^. 
The rotation frequency fl and the magnetic field B are ap- 
plied along the 2-axis. The z-polarized spin current J S (R) is 
induced in the radial direction e r . Here, the spin polarization 
vector is denoted by s. 



with 



and 



eB 

m 



E' = E r 



E c 



(51) 



(52) 



The solution (50) contains two "electric fields" : the spin- 
independent part E r and the spin-dependent one E CT as 
discussed in Sec. IV. The spin-independent part yields 
the conventional Hall effect in the rotating frame. If 
the ends in the radial (longitudinal) direction of Pt film 
attached to the rotating disk (see Fig. 2) are electri- 
cally connected, the Hall voltage is obtained in the az- 
imuthal (transverse) direction while the spin-dependent 
part causes the spin current generation. The z-polarized 
spin current generated by the mechanical rotation can be 
evaluated by 



Js = -eJ dkTrk^/^k)*-] 



(53) 



which leads to the explicit form of the spin current: 

3 S (R) = J r s (R)e r + Jf(R)e (t> , (54) 

with 



j: 



1 + (rw c ) : 



- 7° 



Jt = 



1 + (™ C ) 2 S 



(55a) 
(55b) 



Here, e r (e^,) is the unit vector of the radial (azimuthal) 
direction as indicated in Fig. 2, and is given by 



J s °(i?) = 2nenuj c R. 



(56) 



In the large lo c t limit, the radial spin current, JJ, van- 
ishes and we have 



Jt -> J?. 



(57) 
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This reproduces the ballistic case, Eq. (35). Putting 
lj c t -C 1, the radial spin current bacomes much larger 
than the azimuthal one: 

j; » jf. (58) 

In the following section, we consider the latter case. 



VI. DETECTION METHOD OF THE 
ROTATIONALLY INDUCED SPIN CURRENT 

In this section, we investigate a detection method of 
the spin current induced by rotation using the inverse 
spin Hall effect (ISHE). 54 ' 55 The ISHE converts a spin 
current into an electric voltage by means of the SOI as 



E 



ISHE 



oc J, xs 



(59) 



with the electric field due to the ISHE, Eishe, the spatial 
direction of the spin current J s , and the spin polarization 
vector of the spin current s. In the Pt film attached to a 
rotating disk (see Fig. 3), the electric field Eishe is gen- 
erated in the azimuthal direction since the z-polarized 
spin current 3 S (R) is mainly induced in the radial direc- 
tion when uj c t <C 1. 

Our setup of the Pt film is an open circuit, which means 
that the conventional Hall effect originating from E r in 
Eq. (50) is not created. In the film, the charge degree of 
freedom is quickly frozen due to the charge accumulation 
at the ends of the sample. 

In contrast, owing to the nonconservative nature of 
the spin degree of freedom, spin current is not fully com- 
pensated and can flow steadily even in the open circuit. 
The compensation of spin current occurs within the spin 
diffusion length from the ends. We obtain the spin accu- 
mulation at the sample ends by solving the spin diffusion 
equation with a source term originating from the rota- 
tionally induced spin current. 34 The explicit solution of 
the spin diffusion equation is given in Appendix D . 

The voltage of the ISHE, Vishe, is estimated by 55 



Vishe = ®wp{J r s 



(60) 



with the spin Hall angle O, sample width w, and resis- 
tivity p. The average of the radial spin current (JJ) is 
calculated by 



(Js) 



JldR : 



Ri 



co^TneK(Ri + R%) 

1 + (WcT) 2 



From Eq. (61), we obtain the voltage 

ViSHE OC B 2 n, 



(61) 



(62) 



when uj c t <C 1. This shows that the sign of the voltage 
changes when the rotation axis is reversed whereas it does 
not when the magnetic field is reversed (see Fig. 3). 
In the case of the Pt thin film with w — d — I0mm, 
lu c t = 0.01, 6 = 0.01 and p = 12.8 x 10 _8 fi-m, which is 




B[T] 



-5 5 

Q[kHz] 



FIG. 3. The voltage induced by the inverse spin Hall effect 
in a rotating Pt film is plotted. The Pt film is attached to a 
rotating disk at 50mm from rotation axis. Vishe is propor- 
tional to B 2 and fi when lj c t <C I. (a) Vishe is plotted as a 
function of the magnetic field B at Q. — IkHz. (b) Vishe is 
plotted as a function of the magnetic field Q at B = IT. 



attached to a rotating disk with R ~ 50mm, the average 
of the radial spin current is 

(JJ) » I0 5 A/m 2 (63) 

and the voltage accross the width of the sample is esti- 
mated (see Fig. 3) as 



Vishe ~ l^V. 



VII. LINEAR ACCELERATION 



(64) 



In this section, we discuss spin current generation by 
the linear acceleration in the absence of electromagnetic 
fields. The Dirac Hamiltonian in a linearly and rota- 
tionally accelerating frame without electromagnetic fields 
was derived by Hehl and Ni: 39 

H =/3mc 2 + cot ■ p -\ [(a • r)(p • a) + (p • a) (a • r)] 



2c 

+/3m(a • r) - • (L + S) 



(65) 



with the linear acceleration a and the angular momen- 
tum L = r x p. In the low energy limit of this Hamilto- 
nian up to the order of 1/m 2 , one has the Hamiltonian 
of the Pauli-Schrodinger equation of electron's 2-spinor 
wave function in the accelerating frame: 39 



H 



P 

2m 



4mc 2 



ma ■ r - fi • (L + S) 

a ■ (a X p), 



(66) 



where the rest mass energy and the redshift effect of the 
kinetic energy 39 are neglected. We note that the last 
term, the so-called "inertial spin-orbit interaction" , does 
not contain the mechanical rotation because of the ab- 
sence of the magnetic field. It should be emphasized that 
the inertial effect of the linear acceleration on electon can 
be interpreted as an "effective electric field" . Introducing 
the "electric field" 



E„ 



(m/g)a, 



(67) 
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the inertial SOI can be rewritten as 



A 



(a) 



ffs,a = -tr • (p X qE a 



(68) 



Together with the second term of Eq. (66), ma-r = <7E a - 
r, the inertial effects of the linear acceleration without 
rotation on the large spin-orbit interacting system can 
be analyzed by the same framework as the conventional 
Hamiltonian with SOI: 



2m h 



dU 
dr 



x p 



(69) 



with the potential U = q~E a ■ r. 

The electron velocity in the linearly accelerating frame 
is given by 



where 



P 

r = — 

TO 



v CT . a = ~r<y x E a 



(70) 



(71) 



Let us focus on the inertial effects of linear accelera- 
tion. The anomalous velocity v CT a yields the mechanical 
analogue of the spin Hall effect in a ballistic regime. The 
i— polarized spin current (i = {x, y, z}) generated by the 
linear acceleration is estimated as 



« m r i 2nem\ 
J s = enTr[o-iV ., a J = — - — Sj 



x a. 



(72) 



When the acceleration is induced by the harmonic oscilla- 
tion with the frequency uj a and amplitude u in a;-direction 



uuj a e a e T 



(73) 



the z-polarized spin current is created in y-direction 



2neTO,A 



mole 1 ' 



(74) 



Assuming the Pt film vibrated with ui a — 10GHz and 
u = lOnm, the ac spin current is estimated to be J| ~ 
10 7 A/m 2 (see Fig. 4). It is a future challenge to probe 
the ac spin current in the high frequency mechanical 
resonator. 57 In such a noninertial system, it is straightfor- 
ward to extend the results in the ballistic regime to those 
in the diffusive regime using a well-established framework 
on the spin Hall effect by replacing the usual electric field 
E with the effective one E a . 



VIII. CONCLUSION 

In this paper, we have investigated theoretically the 
generation of spin currents in both rotationally and lin- 
early accelerating systems. The explicit form of the spin- 
dependent inertial force acting on electrons in a rotating 
frame in the presence of electromagnetic fields was de- 
rived from the generally covariant Dirac equation. It has 



J, 



Pt 



„2 iaj 

a — ua>„e " 




coJGHz 



i/[nm] 



FIG. 4. Spin current generation in a linearly accelerating 
frame is shown, (a) When a Pt film is attached to a mechan- 
ical resonator and vibrated in rc-axis, the z-polarized ac spin 
current is created in the y-direction. (b) The amplitude of Jl 
is plotted as a function of o; a at u — lOnm. (c) The amplitude 
is plotted as a function of u at oj a = 10GHz. 



been shown that the force is responsible for the genera- 
tion of spin currents by mechanical rotation in the first 
order of the spin-orbit coupling. The effect of fluctuation 
of the rotation axis on the spin current was discussed 
using the time-dependent part of the force for future ex- 
perimental analysis. 

We have also studied the spin current generation from 
the mechanical oscillation. The spin current can be cre- 
ated in a uniformly oscillating conductor with a large 
spin-orbit coupling because of the inertial spin-orbit in- 
teraction originating from the linear acceleration. We 
provided a concise interpretation of an inertial effect of 
the linear acceleration on electron as an effective electric 
field, which allows us to use the conventional theory of the 
spin Hall effect for describing the spin current generation 
due to the linear acceleration, by the simple substitution 
of the effective electric field for an ordinary one. 

The framework proposed here offers a new route to 
study the inertial effects on electron transport phenom- 
ena, leading to an innovative combination of microelec- 
tromechanical systems (MEMS) and spintronics. 
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Appendix A: Metric, vierbein, and spin connection 
in a rotating frame 

In Sec. II, we use the vierbein representation of the 
Dirac equation. The relation of the matrices between 
curved space-time and flat one is provided by the vier- 
bein (or tetrad) filed ^f a \(x), whose indices /i and (a) 
labels the curved space-time coordinates and the local 
flat space-time coordinates, respectively. The vierbein 
is a local orthonormal base {e^ a -j} a= o,i,2,3 and a kind of 
square root of the metric tensor g^ v (x): 

< a) (aO<%) = <TO*0, (Al) 

with the Lorentz metric rf ,v = diag(— 1, 1, 1, 1). Starting 
with the metric tensor: 

5oo = -1 + m 2 , g 0l = gio = Ui, g tj = <%, (A2) 

its inverse tensor is 

9°° = -1, 9° l = 9'° = u u gv = <J« - u iUj , (A3) 

where Ui = (flxY/c) i and u 2 — ul+u^+u 2 . The vierbein 
in the frame is obtained from Eq. (Al) and (A3): 



e (o) = ~ U *> e \j) = S j- 



(A4) 



We also have the inverse of the vierbein e^\x) — 



e — °0 + V u n e i — °i 



(a) _ m 



(A5) 



The Clifford algebra in the curved space-time can be 

written as 7** (a:) = g^ v (x)ei?\x)^ a with 70 = i/3,7» = 
—ifiai. Thus, we have the algebra in the rotating frame 

7° 0*0 = -7o = -iff, 

Y(x) = Ui7o + 7i = ifiui - if3cti, (A6a) 

where 



/ O 

O -I 



/?=[,. r I • « 



O a 
a O I' 



(A7) 



In the vierbein representation, spin connection T^x) 
is expressed as 

r,(x) = -\w e^g»\d»e^ - i>(f)), (A8) 
where the affine connection is defined by 

r£„ = \9 Xa {dv9^ + d»g av - %9^)- (A9) 
Substituting (A2) into (A9), 

pO T^O T^O j^i pi 

1 00 — 1 iO — 1 ij — 1 jk — u > 



1 00 — 



+ d oUt , r;. = -^i. (Aio) 



From the affine connection, Clifford algebra, and vier- 
bein in the rotating frame, we arrive at the spin(spinor) 
connection: 



r (aO 

Ti(x) = 0. 



4c 



2c 



(Alia) 
(Allb) 



Appendix B: Foldy-Wouthuysen-Tani transformation 

In Sec. II, we derive Pauli-Schrodinger equation in 
a rotating frame by Foldy-Wouthuysen-Tani transforma- 
tions (FWTTs) . The details of the derivation are shown 
below. FWTT is a unitary transformation and a block 
diagonalization method of the Dirac Hamiltonian. First 
of all, we define the even and odd part of the Hamiltonian 
as 



H = /3mc 2 +S + 0, 



(Bl) 



with 



£ = qA a - ft ■ (r x 7T + £), (B2a) 
O = col ■ 7T. (B2b) 

The even part, £, is the block diagonal part of the Hamil- 
tonian whereas the odd part, O, is the block off-diagonal 
part. FWTT is defined by 

H' = UHU^ - UdtU^ (B3) 

with 

( iPcOL-Tz" 



U = exp 



2mc 2 ) P ^ 2toc 2 

A low energy expansion of the Hamiltonian H' into an 
exponential series of 1/m gives a systematic expansion in 
a way that odd parts of H' vanish in any order. Up to 
the order of 1/m 2 , Eq. (B3) is reduced to 

O 2 ' 



(B4) 



H' =13 



mc 



2mc 2 



8m 2 c 4 



0,[0,£] 



Using the relation 

aiCtj = Sij + ieij k a k , 
[7Tj,7r,-] = ihqe ijk B k 

the second term of Eq. (B5) becomes 

po 2 



PaiajiTiiTj _ 



TV 

2m 



qh 
2m' 



a B 



2mc 2 2m 

Next, we focus on [O, [0,£] + ihO]. 

[0,£] +ih6 

= {cat ■ tv, — qA ] + col ■ (—iqh)d t A 

+ [col ■ TV, -ft • (E + r x 7r)] 
= -ificqoi ■ (-VA) + d t A) 

—ihea ■ (O x 7r) + ihea ■ (ft x n 

—ihcq(fl x r) x B 
= ihcOL ■ q~E' 



ihO , 
(B5) 

(B6a) 
(B6b) 

(B7) 



(B8) 
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with 

E' = -VA + d t A + (fi x r) x B. (B9) 
Thus, we have 
1 



O, [0 1 £}+ihO 

1 



{cat ■ 7r, —iKca ■ qE/] 



ifi 



qh 2 



+ 



8m 2 c 2 
qh 



divE' 

a ■ (tt x E' - E' x tt). (BIO) 



From Eqs. (B7) and (BIO), we obtain the Hamiltonian 
i/pR, in Sec. Ill, by neglecting the rest energy term. 



Appendix C: Electromagnetic fields in a rotating 
frame 



In this Appendix, we give a brief review of a general 
relativistic transformation of electromagnetic fields. 

Electromagnetic fields in the rotating frame are related 
to those in the rest frame as 



E' = E + (ft x r) x B, 
B' = B, 



(Cla) 
(Clb) 



when the rotation velocity is much less than the speed of 
light, \fl x r| -C c, where E' and B' are electromagnetic 
fields in the rotating frame. 

Equations (CI) are not Lorentz transformation in spe- 
cial relativity but a general coordinate transformation in 
general relativity. The Lorentz transformations of the 
electromagnetic fields are written as 43 



E"/c = 7(E/c + /3 x B) 
B" = 7 (B - 3 x E/c) 



7 



7 + 1 

7 2 
7 + 1 



(J3 ■ E/c)/3 (C2a) 



(6-B)8, (C2b) 



with 



7 = \ = ^ 



(C3) 



Here E" and B" are the electromagnetic fields in the iner- 
tial frame which has a uniform velocity Vo relative to the 
rest frame. Replacing Vo with v(x) = ft x r in Eqs. (C2) 
can never reproduce the correct transformation between 
the rotating frame and the rest frame (CI). The special 
relativistic transformations (C2) have an apparent sym- 
metry for E/c and B whereas the relations (CI) do not. 
Such an asymmetry in (CI) originates from the space- 
time asymmetry of a general coordinate transformation 
which relates physical quantities in a rotating frame to 
those in a rest frame. 



When the rotation axis is parallel to the z-axis, the 
transformation between the rest frame and the rotating 
frame is expressed as 



with 



( 1 \ 

cos fit sin fit 

—sin fit cos fit 

\0 1/ 



(C4) 



(C5) 



where rotating coordinates carry the prime. Equation 
(C4) leads to a general coordinate transformation 42,43 



R = 



1 

dx a | fly/c cos fit sin fit 

dx'i 3 ~ | -flx/c sin fit cos fit 

1 



(C6) 



According to the principle of general covariance, electro- 
magnetic fields are components of a second rank tensor 
in general coordinate transformations. Electromagnetic 
tensors in the rest frame F and those in the rotating 
frame F' are related as: 



F' = LR T FRL T , 
in which T denotes the transpose matrix, and 



/ -E x /c 
E x /c 
E y /c -B z 
V E z /c By 



-Ey/C 

B z 




-Ez/c \ 

-By 

B x 




(C7) 



(C8) 



Equations (C7) and (C8) lead to 



K 


— E x 


+ xQB z 


(C9a) 


K 


= E v 


+ yQB z 


(C9b) 


K 


= E Z 


- n(xB x +yB y ) 


(C9c) 


K 


= B x 




(C9d) 


K 


= B y 




(C9e) 


B' z 


= B Z . 




(C9f) 



Thus, we obtain the relations (CI) in Sec. IIIB. The 
relations shown in this section hold for \ft x r <C c. 
Generalized relations for electromagnetic fields in more 
rapidly rotating frame should include a Lorentz fac- 
tor 7 = l/\/l — (ft x r/c) 2 . 42 We omit discussion on 
Maxwell's equations in a rotating frame. The subject 
is originally studied by means of generally covariant 
Maxwell's equations 44 . 



Appendix D: Solution of the spin diffusion equation 
in a rotating frame 

In this Appendix, we discuss the spin diffusion equa- 
tion, to clarify the spin-dependent transport property of 
spin current in a rotating frame. In the non-equilibrium 
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steady state, spin current in the rotating frame consists 
of the diffusive current jf m and the drift currnt Jl ot : 



J, 



tot 



4 m + ji 



with 



^ = _ V(5 

ep 



(Dl) 



(D2) 



Here, dp is the spin accumulation. The drift part, which 
can be regarded as a spin source term, is given by the 
mechanically induced spin current (54) in Sec. VI. The 
spin diffusion equation with the source term is written as 



W 2 Sp 



1 

a! 



Sp — epdivJ s 



(D3) 



with the spin diffusion length A s . When both rotation 
axis and magnetic field are parallel to the z-axis, J™* is 
mainly induced in the radial direction in the Pt film as 
discussed in the text. Thus, we solve this equation in 
the cylindrical coordinate. The radial component of the 
equation is reduced to 



1 d ( dSp\ Sp ep d 



r dr 



dr J A 2 , r dr \ 1 + (tlo c ) 2 



,- jggg£).(D4) 



Using parameters 



a = — , b = 2epA, 

A s 



A 



2nenTUj 2 
1 + (™ c ) 2 



(D5) 



the analytical solution can be written as 



5p(r) = ^2 + c i Jo(ar) + c 2 K (ar), (D6) 



where J n , K n are the modified Bessel function of the first 
and second kind, and the constants c\ and c 2 are to be 
determined by the boundary condition: 



Jf t (i? l ) = 0, (i = l,2) 



(D7) 



at the ends of the Pt, r = R\, R 2 . In this condition, we 
have 

ci = c [i? 2 Ji(aRi) + R 1 J 1 (aR 2 )] (D8a) 
c 2 = c [-R 2 K 1 (aR 1 ) + R 1 K 1 {aR 2 )} (D8b) 

with 



c = 



a[J 1 (aR 2 )K 1 (aR 1 ) - J 1 (aR 1 )K 1 (aR 2 )] 



• (D9) 



Putting d = lOOnm, R x = 50mm, Q, = 1kHz, B = IT, 
A s = 14nm, and p = 12.8 x 10" 8 O-m for Pt, 56 we obtain 
the spin accumulation near the ends, Sp(Ri) « 0.05nV 
and Sfj,(R 2 ) w -0.05nV. 
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